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S^ . The self-similar cluster fluctuations of directed bond percolation at the 

percolation threshold are studied using techniques borrowed from intermit- 
tency-related analysis in multi-particle production. Numerical simulations 
based on the factorial moments for large 1 + 1-dimensional lattices allow 
to handle statistical and boundary effects and show the existence of weak 
but definite intermittency patterns. The extracted fractal dimensions are 
in agreement with scaling arguments leading to a new relation linking the 
intermittency indices to the critical exponents and the fractal dimension of 
directed percolation clusters. 
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1 Introduction 

The transition between integrable and chaotic behaviour has attracted inter- 
est for quite a long time. In particular, in the context of hydrodynamics, the 
spatio-temporal intermittent behaviour has been studied [|J. On the other 
hand, similar structures have shown up in the context of particle fragmenta- 
tion at high energies, as first exemplified by the JACEE cosmic-ray event and 
later on by many studies at particle accelerators @. Intermittency studies 
have also been performed in Statistical Mechanics, on systems like the 2D 
Ising model at the second-order transition || f|. In this work, we consider 
directed bond percolation as a simple toy model for intermittent behaviour. 
As we shall argue below, it might provide a convenient model approach to 
describe some basic properties of intermittent systems. 

Our interest in studying intermittency was triggered by the observation 
that directed percolation can be used as a simple non-perturbative model 
for particle fragmentation as described by the Schwinger mechanism [|5], 0]. 
Already before that, it had been proposed that directed percolation might 
be useful in describing the transition between laminar and turbulent flow, 
occuring for example in hydrodynamics JT|. We shall describe these heuris- 
tic arguments in section 2. For the quantitative description of intermittent 
behaviour, a practical tool are the factorial moments 

p ._ = <n(n-l)---(n-q+l)> 
q ' <n>i l ' ' 

where n is the number of particles in a certain (rapidity) interval and <> 
denotes the average over a given sample. We shall recall some of their prop- 
erties in section 2. Intermittent behaviour means that the F q obey a power 
law as a function of the binning size S 

F q {6)~6-t"> (1.2) 

which defines the intermittency indices <p q . The verification of this relation 
and the eventual calculation of these numbers will be the first objective of 
this paper, to be described in section 3. From the 4> q , one obtains the Renyi 
codimension V q 

T> q = <j> q /{q-l) (1.3) 

If the fluctuations of the model have a fractal nature, the T> q will be q- 
independent, while for a multifractal, the T> q retain a dependence on q. In 



the case of a simple fractal behaviour, when T> q = T> for all q, we shall obtain 
a relationship with the fractal codimension D of the directed percolation 
cluster and critical exponents (see section 4) 

V = (3/v L -D>0 (1.4) 

This result allows to reconsider the recent comparison between numerical 
calculations of fractal cluster dimensions with theoretical ideas ||. It also 
casts a bridge with the Reggeon Field Theory, which furnishes a continuum 
field theory description of directed percolation || . The validity of the result 
(1.4) only depends on simple scaling arguments and should therefore not be 
restricted to critical directed percolation systems. Section 5 presents our 
conclusions. 



2 Directed percolation and intermittency 

2.1 Some heuristics 

We recall the definition of directed percolation (see e.g. [|Kj). Consider a 
lattice and single out a preferred direction. This preferred direction will be 
called "time", the other ones are referred to as "space". For the purposes of 
this paper, we shall limit ourselves exclusively to the case of one temporal 
and one spatial dimension. The sites of the lattice can be either filled (1) 
with probability p or empty (0) with probability 1 — p. A bond between 
sites is present with probability p and absent with probability 1 — p. The 
occupation v^t of the sites is defined in terms of the conditional probabilities 
P(vi,t+i\vi-i,uVi+i,t), see fig. 1, namely, 

P(1|0,0) = 0, P(l\l,l)=p(2-p)p , P(l|0,l) = P(l|l,0)=pp (2.1) 

together with P(0\v,v') = 1 — P(l\v,v'). For p = 1, one has directed bond 
percolation and for p = 1 one has directed site percolation. The probability 
V for having an infinite directed cluster of connected sites scales with the 
distance from the percolation threshold p c as V ~ (p — p c ) 2 ^ and has the role 
of an order parameter. The scaling of the temporal and spatial correlation 
lengths £|| and £j_ is described by exponents v\\ and v>±, viz. 

£||~(P-Pc)- V|1 , iL~(p-Pc)-" X 



Numerical values are ill 



p c = 0.705489(4) square site 

p c = 0.644701(1) square bond (2.2) 



and for the exponents JTTJ] 

P = 0.277(2) , u\\ = 1.7339(3) , v± = 1.0969(3) (2.3) 

From these, all other exponents can be found by scaling relations. For a 
test of universality for the exponents and the two-point correlation scaling 
function, see [^ . 

The problem of directed percolation can be rewritten in terms of a field 
theory. As it is well known, directed bond percolation is in the same uni- 
versality class as Reggeon field theory ||, at least close to the upper critical 
dimension D* = 5. Nevertheless, this universality was also confirmed numer- 
ically for D = 2 and D = 3, see refs. |L3], [14], |15| . Consequently, results of a 
directed percolation study also carry an implication about that theory. We 
will discuss that point later on. 

Some time ago, Pomeau [fj has proposed that spatial-temporal intermit- 
tency were in the same universality class as directed percolation. According 
to him, the transition to the turbulent state in a liquid is similar to the one 
induced by the intermittent transition in oscillations and that the interaction 
between neighboring oscillators might be seen as some sort of local contam- 
ination process. It appears then plausible that an oscillator in a turbulent 
state may either relax spontaneously towards its quiescent state or else may 
contaminate its neighbors. As an example of this, one might consider the 
observation that the laminar convection rolls in Taylor-Couette flow do not 
break up spontaneously into a turbulent state but get only infected via a tur- 
bulent neighbor. This is nothing else than the familar definition of directed 
percolation. Although this argument is purely qualitative and does not yield 
a value for the percolation probability p, it does carry some heuristic value, 
even if the simple directed percolation is probably not sufficient to describe 
intermittency occuring in realistic liquids Jl6| . 



To complement this, we now argue that directed percolation might be 
useful to study some aspects of particle fragmentation. For simplicity, con- 
sider the example of e + e~ annihilation producing a quark-antiquark pair. It 
is generally believed that for a confining theory like QCD the color lines of 



force are essentially concentrated in a narrow tube connecting q with q, act- 
ing like a string with a constant tension. By the Schwinger mechanism, new 
qq pairs are created out of the vacuum by the color force with a probability P 
expressible in terms of known quantities || |6| . One thus obtains a resulting 
color field which will strongly fluctuate as a function of the relative rapidity 
between the particles. 

We now make the link from directed percolation to this model. Since 
the color flux tube is essentially one-dimensional, we take a two-dimensional 
lattice. The direction perpendicular ("space") to the imposed propagation 
("time") direction will be taken to correspond to the rapidity interval. On 
each site of the spatial lattice we put if there is no color field and 1 if there 
is one. The quarks or antiquarks generated are placed on the dual spatial 
lattice. The creation of a new qq pair generates a supplementary color field 
which exactly compensates the one which produced the new pair. In the par- 
ticle language, this means that some local field energy has been transformed 
in particle creation with probability P. In the percolation language, this 
corresponds to having a site occupied by at time t after its two neighbors 
at time t — 1 have been occupied by one. In the frame of directed bond 
percolation, see (2.1), we have P = 1 — P(l|l, 1) = (1 —p) 2 - Again, this very 
simple model should be taken as a guide but not as a realistic description of 
high-energy processes. 

The breaking of the color field is stochastic and, due to the causality 
nature of directed percolation, new particles can only be created where the 
field persits, it generates a cascade. Since randomness and cascading are 
just the two features known to give intermittency in particle production 
[§], and because of the connection with turbulent behaviour [|l], we expect 
to find intermittency in the context of directed percolation. This means the 
occurence of short-range self-similar fluctuations which dominate the average 
characteristics of a multiparticle system. 

This rough model description motivates our simulations of 1 + 1 dimen- 
sional directed percolation for particle physics. However it is worth noting 
that 2 + 1-dimensional directed percolation has another more theoretical re- 
lation with multi-particle collisions, thanks to its continuum limit in terms 
of the Reggeon Field Theory. In this context, the "time" dimension corre- 
sponds to the rapidity of produced particles in a very energetic collision, and 
the two-dimensional "space" is the transverse impact parameter. So, some 
of our results can find a prolongation in a 2+ 1-dimensional study of direct 



relevance to the theory of particles. 

2.2 Factorial moments 

The factorial moment analysis |7j is useful to avoid as much as possible 
a "contamination by statistical noise" of the dynamical fluctuations one is 
looking for. The general idea is that a large part of the non-dynamical 
statistical fluctuations is due to the small number of particules (statistical 
realizations), and thus can be assimilated to a "noise" with simple stochastic 
characteristics. In the absence of other information, the "noise" is assumed to 
be Poissonian and thus can be de-convoluted from the observed distribution 
of particles by the factorial moments' trick. Of course, this assumption looks 
strong, but many model tests [[71 have shown its validity in multi-particle 
physics. On the other hand, statistical systems like directed percolation are 
useful to check the "noise" assumption in cases where it is embedded in the 
dynamics, and allow to deepen the statistical analysis. 

In this specific case, however, one disposes of an information on the total 
number of sites under study in a given phase-space box, and thus, adding 
this constraint to the "noise" definition, one is led to consider the Bernoulli 
"noise". More specifically for lattice systems, one writes : 



F q 



<p q > 
< p >i 



q = <n(n-l)...(n-q + l)> 
P ~ <N{N-l)---{N-q-l)> [ ] 

where N stands for the number of sites per bin, p is the density of occupied 
sites and the average <> is meant over the total number of bins. Notice that 
the choosen normalization for factorial moments has been designed for the 
Bernoulli noise J7J expected in lattice simulations [Hj] and checked in detailed 



calculations on the Ising model |L8|. If the factorial analysis works well, the 
density p and its fluctuations are free from statistical perturbations and the 
intermittency parameters 4> q and V q (see (1.2), (1.3)) can be estimated. 

Intermittency studies for statistical systems using factorial moments (2.4) 
are particularly interesting since they offer an unique opportunity to discuss a 
scaling behaviour of fluctuations in a well-defined dynamical scheme. There is 
no need of a small coupling hypothesis nor of an assumption on the statistical 



noise. Indeed, the statistical fluctuations from site to site are, together with 
the dynamical ones, encoded in the generic rules (2.1). 

3 Numerical simulation 

We now turn to describe our simulation procedure and the numerical results 
for the intermittency indices. We use directed bond percolation on a square 
lattice. Since an intermittent behaviour as the one in eq. (1.2) can only occur 
at criticality, we shall always work at the critical point p = p c , see (2.2). As 
the initial state, we take occupied sites (1) on all spatial sites at time t — 0. 
Both periodic and free boundary conditions will be considered, see below. 

The factorial moments can be calculated for a two-dimensional or a one- 
dimensional binning. In the latter case, we can either choose a spatial (verti- 
cal) or a temporal (horizontal) one. For the calculation, one can either take 
subdivisions of the entire lattice with N spatial sites into account or only 
apply the procedure of repeated subdivions to a certain subcell of L < N 
spatial sites. Both of these procedures were tried for the case of periodic 
boundary conditions, while for free boundary conditions, only the second 
one was used. The moments were obtained calculating by generating L time 
steps of the evolution of the directed percolation cluster. The resulting square 
of L x L sites was then subdivided into squares of sizes L/2 b x L/2 b with 
b = 0, 1, 2, ... to obtain the two-dimensional factorial moments F^ 2D \ To ob- 
tain the one-dimensional moments F^ and F^ v \ in the same square for the 
fixed time (for F^) or for fixed position (for F^ v ">) subcells of length L/2 b 
were formed, see fig. 2. We take an increment of five "layers" (intermediate 
steps in time or space, respectively) between two consecutive calculations 
for the ID factorial moments with a given square of size L x L. To re- 
duce the statistical correlations between measurements, we generated 2000 
"thermalisation" time evolutions steps before beginning the calculation of 
the factorial moments and generate an additional 2000 time steps between 
any two squares to be considered for measurements. Also, if the directed 
percolation cluster happened to die out (this only occured for fairly small 
lattices), we reinitialized with all spatial sites occupied and 2000 thermalisa- 
tion steps. Since correlation function studies have shown that scaling sets in 



after about 20 thermalization steps (see [12 1), this thermalization should be 



sufficient. Almost the entire CPU time was spent in calculating the factorial 



moments. 

We now discuss the necessary spatial lattice size in order to become inde- 
pendent of finite-size effects. To see this, we examine the calculated factorial 
moments F^S) as function of the binning size S. In fig. 3-a, we display 
the Renyi codimension X> 4 as calculated for boxes of linear size L/2 b for the 
2D factorial moment F 4 . Similarly, fig. 3-b shows the same quantity as 
obtained from F 4 and fig. 3-c those obtained from F 4 . 

We first observe that for all three cases the apparent Renyi codimen- 
sions become independent of the lattice size N, the size L of the largest box 
considered and of the boundary conditions if N is larger than 4096. These 
sizes are somewhat larger than the largest ones considered in intermittency 
model calculations so far (see e.g. [[18~1). We can conclude that the data for 
N = 8192 will describe accurately the behaviour of the N — > oo limit and 
that finite-size effects are small for this size. Secondly, we see that, apart 
from the largest blocks (where b is small) and the smallest one (where the 
cells just contains two sites), the data show a large and wide plateau and 
become b- independent. The value of the plateau will be taken as our esti- 
mate for the Renyi codimensions. These observations hold for all the three 
types of factorial moments considered. We have drawn similar plots for the 
cases q = 2 and q = 6 and find the qualitative behaviour to be the same, 
although the convergence towards the N — > oo limit is somewhat slower for 
higher rank. We can conclude that our directed bond percolation model at 
p = p c does indeed show intermittent behaviour. For a precise calculation of 
the intermittency indices, we have seen that lattices iV = 8192 spatial sites 
are large enough. We have also studied the effect of varying p away from 
p c . Already for variations of order 5p ~ 10~ 2 , the T> q are no longer almost 
independent of b, but are rather seen to rise steeply with increasing values 
of b. 

To investigate the q- dependence, we take N = 8192, L = 2048 with 
periodic boundary conditions. We generated 500 squares of size L x L for 
the calculation of the F q , with 2000 intermediate time steps in between each 
square. In fig. 4-a, we give the curves for the Renyi codimension T> q 2D \ for 
q = 2,3, 4, 5, 6 as calculated from FJ? D >. Similarly, fig. 4-b shows the results 
obtained from Ff> H ^ and fig. 4-c those from FJ- V \ Already a first glance at 
the figures shows that the 2D and the horizontal (H) results are very similar 
to each other while the vertical (V) ones are quite distinct. In both figs. 4-a 



and 4-b, we see that for a large range of binning boxes (as given by b), the 
apparent Renyi codimension T> q , is in fact q- independent. This means that 
the 2D and the horizontal intermittency behaviour is in fact described by a 
simple fractal structure. We read off the final Renyi codimensions 

p( 2D ) = 0.023 ± 0.004 (3.1) 

for the 2D factorial moments and 

V {H) = 0.021 ±0.004 (3.2) 

for the horizontal factorial moments. In section 4, we shall return to an 
interpretation of these results. 

On the other hand, the data obtained from the vertical factorial moments 
show quite a distinct behaviour, see fig. 4-c. It is not possible to define 
anything like a single Renyi codimension. To our best knowlegde, this is the 
first example of a system displaying such a strong anisotropy in the factorial 
moments. We note that for the maximal values of b, the effective T) q are 
close to the values obtained from the 2D or the horizontal case. This is to 
be expected since for these very small boxes the effects of the anisotropy of 
the directed percolation clusters have not yet become apparent. 

4 Scaling relations and intermittency indices 

The features revealed by our numerical simulations and the estimates (3.1) 
and (3.2) are quite unexpected and thus require some theoretical interpreta- 
tion. In particular, the obtained intermittency dimensions are much smaller 
(while non-zero) than the exponent ratios (3/v\\, P/v±, (see (2.3)), contrary 
to some expectation in non-directed percolation @]. Moreover, the similarity 
between p( 2D ) and T>^ has to be explained. We shall now do so by using a 
simple scaling argument, eventually arriving at eq. (1.4). 
Let us first consider the two-dimensional second moment: 

Ff D) :=< n(n - 1) > / < n > 2 

where one keeps the simple normalization (1.1) which is sufficient for the 
theoretical discussion in the continuum limit. As already noticed in the 
literature [EL H], the intermittency behaviour (1.2) is to be related to an 



effective singularity in the relevant correlation function. In the case of the 
2-D moment, one may write: 

< n(n — 1) > oc / dtidt 2 dxidx 2 G^ 2 \xi,x 2 ;ti,t 2 ) (4.1) 

where G^ is the two-point correlation function and the integral is over the 
phase-space box of unit length 5, see fig. 2. This corresponds to the proba- 
bility of a connecting path of occupied sites between the space-time points 
(xi,ti) and {x 2 , £2)- Note an implicit averaging over phase-space boxes which 
is allowed by translational invariance, well verified in our simulations. It is 
also important to notice that, for practical simulation purpose, the binsize S 
has to be considered very small, but large with respect to the bond size. 

The scaling behaviour of G^ near the percolation transition point p c can 
be written: 

28 

G {2) ~ (5p) 2(i V[(6p) v± 6; (<*?)"" r] ~ 9~~ ■ cste (4.2) 

where 5p = p — p c , 9 = \x± — x 2 \ and t — \ti — 1%\ and \I/ is a scaling function. 
The first part of relation (4.2) comes from the critical behaviour at the 
percolation transition while the second part is the dominant <5p-independent 
behaviour after rescaling. Note that the known strong anisotropy in time of 
directed percolation provides a hint for the dominance of the 9 singularity 
over the phase-space box. This typical feature of directed percolation is 
well demonstrated by the narrowness of percolation clusters in the "time" 
direction, see ref. 0. Inserting the behaviour (4.2) in the estimate (4.1) at 
small S one finds: 

< n(n - 1) >~ 5 4 - 2 ^^ (4.3) 

the exponent 4 being simply interpreted as the naive integration dimension 
in (4.1). 

The key remark of our analysis now coming is related to the normalization 
denominator < n > 2 of F 2 . An oversimplified evaluation of this denomina- 
tor would have led to some number proportional to the phase-space, namely 
5 4 , and thus to an intermittency exponent /3/v±. This is much too large 
as compared with the numerical value (3.1). Rather, one has to take into 
account the density of occupied sites in the box [6], which is not uniform. 
This density is better described by a fractal set whose dimension D has been 



recently numerically estimated from the formation of very long percolation 
clusters 0. Using this information, we write: 

< n >~ 5 D = 5 2 ~ D (4.4) 

introducing the fractal dimension of occupied sites and its codimension D. 
We thus get the general relation: 

if D) ~ ( r 2 ) (/3/ ^- 5) (4.5) 

Comparing with (1.3), we note that since the boxes considered are two- 
dimensional, we have F^ D > ~ (6~ 2 )^ q . In terms of the Renyi codimensions 
p( 2D ) of a (fractal) intermittent process, we recover indeed the relation (1.4). 
The factorial moment satisfies 



i(2D) . 



< n(n - 1) > < p > 



<n> 2 < p > 2 

which implies in turn that the Renyi codimension T> is non-negative. From 
our numerical estimate (3.1), we see that the fractal cluster codimension D 
and the exponent ratio j3 /v±_ are very close, up to a small difference of definite 
sign. 

It is of interest to compare this result with earlier determinations of the 
fractal dimension of directed percolation clusters. A careful study of this 
question was presented in . Three possibilities to calculate D were consid- 
ered. The first one simply assumes that the two-dimensional fractal can be 
considered as a direct product. The second possibility relies on a determin- 
istic growth argument. Finally, the third possibility proceeds along a scaling 
argument taking into account the presence of the two length scales given by 



£|| and £j_ |19| . These arguments lead to | 



P (— + i] ~ 0.412 direct product 
i 2(3/ (u\\ + v±\ ~ 0.196 deterministic growth \ > 



(3/v± ~ 0.253 scaling [19] 

Via (1.4), this would lead to "D (2D) ~ -0.16,0.06 and 0, respectively. A 
precise numerical calculation using box-counting arguments yields [§] 

D = 1.765 ± 0.005 (4.8) 



and one notes a discrepancy with all of the relations ( |4.7| ). 

However, (4.8) would yield via (1.4) a value £>( 2D ) = 0.018(5) in agreement 
with our estimate (3.1). In other terms, one may interpret the small but 
non-zero difference between the cluster and correlation fractal dimensions 
as due to a weak intermittency phenomenon in 1 + 1-dimensional directed 
percolation. 

(H) 

The same relation (1.4) holds for the "space" moments F 2 , but it re- 
quires a more refined way of derivation, since the 2-point correlation function 
is zero between 2 sites with no "time" separation, see fig. 2. In this particular 
case, one has to consider a 3-point correlation function G^ linking the two 
given sites (xx,ti) and (x 2 ,t 2 ), starting from an initial (x ,t ). One would 
assume the scaling behaviour 

G (3) ~ (5p) 3 ^ [(5p)^6 tJ ; (Sp^nj} ~ 6 12 ^ ■ cste (4.9) 

where 5p = p — p c , 9ij = \xi — Xj\ and r^- = \U — tj\ and i,j — 0, 1, 2;i ^ j. 
Integrating over the initial conditions at (x ,t ), and the phase-space box of 
linear dimension 5 for (xi,ti) and (x 2 ,ti) (note that we are at equal time), 
one gets: 

< n{n - 1) > ~ 5 d ± +i-p/v ± 

<n> ~ 5 D± 
F W „ 5 -Wu ± -D ± ) (41Q) 

where D± is the "horizontal" fractal dimension. One finds again the same 
type of relation (1.4), which thus appears quite general. 

We also note that our numerical estimates (3.1) and (3.2) suggest that the 
fractal codimensions Z)( 2D ) ~ D±, which is in agreement with the literature 



IE 



A simple relation as (1.4) would have to be modified in the case of the 
"time" dimension, due to the multifractal character manifested by our numer- 
ical simulation. It is to be noticed that such a feature is in general associated 



with non-equilibrium cascading processes |2(J. It is thus to our knowledge, 
the first time that multifractal fluctuations appear in a percolation prob- 
lem, in probable relation with the oriented "time" direction characteristic of 
directed percolation. 



5 Conclusions and outlook 

Let us summarize our results: 

1. On the basis of a simulation of directed bond percolation in 1-" space" -1- 
"time" dimensions and after checking against statistical and boundary 
effects the stability of our results, we have found that intermittency 
patterns of fluctuations show up in the distribution of occupied sites 
clusters. The fluctuations along the "space" dimension and in two 
dimensions are dominated by weak, fractal, structures of Renyi codi- 
mension V = 0.022 ± 0.004, (the mean value of (3.1) and (3.2)) while 
"time"-like fluctuations appear with a more complicated structure of 
multi-fractal type. 

2. These numerical results can be interpreted as resulting from the scaling 
behaviour of the two- and three-point correlation function in the con- 
tinuum limit at the bond percolation critical point, both corrected by 
the fractal dimension of clusters. The overall result is summarized in 
formula (1.4), which appears well verified for the two recognized frac- 
tal cases, "spatial" and two-dimensional. These two cases were seen to 
have the same intermittency indices. Our result clarifies a discussion 
on the calculation of fractal dimensions of directed percolation in the 
literature ||. 

It is quite interesting to note that, in its generalized form (1.4), the rela- 
tion expressing the intermittency indices (in the fractal case) as a function 
of the correlation indices and the fractal cluster dimension can be applied 
to other problems. Let us mention the case of percolation at a second-order 
phase transition ||, which is also fractal, according to general scaling argu- 
ments p|. An interesting example is the two-dimensional Ising spin system, 
where the relation (1.4) leads to the following prediction: 

f) = j3/u - D = 7/96 (5.1) 

where the exponent ratio fi/u = 1/8 is well known and the value of the fractal 

codimension D = 5/96 can be obtained by conformal invariance techniques 

21|]. The result (5.1), which is in agreement with numerical simulations [|1^ 



and finite-size scaling evaluations [22 , clarifies a set of contradicting issues 



on the subject. While different from the value f> = (3/u, predicted in ref. 

|], the value (5.1) is different from the fractal dimension D of clusters as 

well. Indeed, the system admits a set of fractal fluctuations around a fractal 

object, which here was the average percolation cluster. We also remark that 



the value found in |22[ is not accidental, as conjectured by the authors, but 
a prediction valid at the continuum limit. 

Finally, the access to so refined quantities as fluctuation moments allowed 
by numerical calculations of directed percolation is an incentive to perform 
the same analysis in 2+1 dimensions, with the bonus of a direct application 
to particle physics through the Reggeon Field model. Note that the already 
obtained result in 1+1 dimension, applied to the Swinger-type model |5|, §], 
confirm the existence of weak intermittent fluctuations, seen in the model 
simulations 01. 
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Figure captions 



Figure 1: The evolution as described by directed bond percolation. The 
result of a simulation on a small lattice with 32 " spatial" sites and periodic 
boundary conditions is shown. The triangles show some typical evolution 
cells whose conditional probabilities are recalled below the simulation sample. 
Note that, by convention, white dots are for empty sites and stars for occupied 
ones. 

Figure 2: Binning for the definition of the factorial moments. Typical binning 
boxes are represented on the simulated patterns of sites (same series as Fig.l): 
Two-dimensional, Horizontal and Vertical examples are shown, together with 
typical connected paths which contribute to the singular behaviour at the 
continuum limit (see text). 

Figure 3: Effective Renyi codimension T>2 as obtained from the factorial 
moment (a) .FJ , (b) F 4 , (c)-F 4 as function of the binning size L/2 b . The 
curves corresponds to the following lattices Q : N = 8192, L = 2048 periodic 
boundary conditions (full curve), O : N = 8192, L = 2048 free boundary 
conditions (dotted), O : N — 4096, L = 1024 free boundary conditions 
(dashed) and □ : N = 4096, L = 4096 periodic boundary conditions (dash- 
dotted). Notice the quite larger vertical coordinate scale in fig.(c). 

Figure 4: Renyi codimensions T> q = 4> q /{q — 1) as obtained for N = 8192 
sites, a maximal box size of L = 2048 and periodic boundary conditions 
as obtained from the factorial moments (a) F^ D \ (b) F^ H \ (c) F^ v \ The 
symbols correspond to the following values ofqQq = 2,Qq = 3,Oq = 4, 
□ q — 5 and □ q = 6. Notice the quite larger vertical coordinate scale in 
fig.(c). 
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